NASA Technical Memorandum 101353 


Infrared Surface Temperature 
Measurements for the 
Surface Tension Driven 
Convection Experiment 


Alexander D. Pline 
Lewis Research Center 
Cleveland, Ohio 


( 1C 1353) IHFElflED SOSFBCE 

IfHPEiAIOEE EiASCJEBEBTS FOB £ DEFACE 

lEBSICB DBIVEI CCIViCTICK EIFFEI^M 

lesetft 0*1,., Aug. 

1588 (»A£l) 25 f csCI U8 


G3/35 


M89-2 1224 


Onclas 

0204477 


March 1989 




NASA Technical Memorandum 101353 


Infrared Surface Temperature 
Measurements for the 
Surface Tension Driven 
Convection Experiment 


Alexander D. Pline 


March 1989 


NASA 

National Aeronautics and 
Space Administration 


Lewis Research Center 

Cleveland, Ohio 44135 


Trade names or manufacturers' names are used in this report for identification 
only. This usage does not constitute an official endorsement, either expressed or 
implied, by the National Aeronautics and Space Administration. 



Contents 


Page 

Summary 1 

Introduction 1 

Symbols 2 

Surface Tension Driven Convection Experiment 2 

Basic Description 2 

Science Requirements 4 

Measurement of Absorption Coefficients of Silicone Oil 4 

Absorption Coefficients 4 

Apparatus and Procedure 5 

Results and Discussion of Absorption Coefficient Measurements 6 

Surface Temperature Measurements 8 

Apparatus 8 

Procedure 12 

Results and Discussion 13 

Effective Emissivity of Silicone Oil 13 

Thermocouple Measurements and Calculated Surface Temperature Distribution 13 

Constant Temperature Experiment 14 

Constant Flux Experiment 16 

Concluding Remarks 18 

Acknowledgments 18 

Appendix— Inframetrics Thermal Imaging System Specifications 19 

References 24 


Summary 

In support of the Surface Tension Driven Convection Experi- 
ment (STDCE), a planned space transportation system (STS) 
flight experiment, a commercially available infrared thermal 
imaging system is used to quantify the imposed thermal 
signature along the free surface. The system was tested and 
calibrated for the STDCE with ground-based equivalents of 
the STDCE hardware. Before using the system, consideration 
was given to the radiation characteristics of the target (silicone 
oil). Absorption coefficients were calculated to understand the 
“surface” depth as seen by the imager and the penetration 
depth of the surface heater (C0 2 laser). The performance and 
operational specifications for the imager and image processing 
system are described in detail to provide an understanding of 
the equipment. Measurements made with the system were 
compared with thermocouple measurements and a calculated 
surface temperature distribution. This comparison showed that 
in certain regions the IR imager measurements were within 
± 5 percent of the overall temperature difference across the 
free surface. In other regions the measurements were within 
±10 percent of the overall temperature gradient across the 
free surface. The effective emissivity of silicone oil for these 
experimental conditions was also determined. Measurement 
errors and their possible solutions are discussed. 

Introduction 

Noncontact surface temperature measurement methods have 
a distinct advantage over intrusive temperature measurement 
methods; it may be undesirable, or in some instances 
impossible, to contact the surface whose temperature is to be 
measured. Examples of such situations include dynamic targets 
at substantial distances (flying aircraft, remote earth sensing 
applications), large targets (buildings), harsh environments 
(applications in industrial areas), and processes which are 
sensitive to disturbances (fluid flow in microgravity). 

The Surface Tension Driven Convection Experiment (STDCE) 
is a planned space transportation system (STS) or space shuttle 
microgravity flight experiment to study thermocapillary flows 
in a reduced-gravity environment. Thermocapillary flow is 
driven by a thermally induced surface tension variation along 
a liquid free surface. This surface tension gradient, established 
by heating the surface, is a surface-tractive force that creates 
a flow parallel to the liquid free surface from regions of low 


surface tension to high surface tension. Although natural 
convection generally predominates in a normal earth environ- 
ment, thermocapillary flows in reduced gravity may become 
significant (ref. 1); these flows could, under certain conditions, 
detrimentally impact planned commercial space activities in 
which a free surface exists (e.g., containerless processing). 
The STDCE is a flight experiment to quantitatively investigate 
thermocapillary flows in reduced gravity. 

The theoretical development for these types of flows as well 
as the influence of configuration and parametric changes on 
the flows are discussed in detail in references 2 and 3, 
respectively. From these analyses the thermal signature along 
the free surface emerges as a very important parameter, 
because it is the driving force for these flows. Therefore, 
quantification of the thermal signature is essential for both 
ground-based research and the STS flight experiment. 
Temperature measurement by a contact method is sensitive 
to free-surface position fluctuations caused by motion of the 
free surface. Such measurements may result in inaccurate data, 
especially when the thermal boundary layer along the free 
surface is thin, as in the case of the STDCE. In the micro- 
gravity environment of the STS a liquid/gas free surface is 
not stationary, because of small accelerations aboard the 
vehicle; therefore, a noncontact temperature measurement 
system is sought for both the STDCE ground-based research 
and the flight experiment. 

In recent years, developments in photo-optical instru- 
mentation technology have made high resolution, noncontact, 
infrared temperature measurements possible for an experi- 
menter who does not have the time or resources to develop 
such a system. Even though a basic system is commercially 
available, it must be tested with, and calibrated for, a particular 
application before it can be used with confidence. Consid- 
eration must be given to the radiation characteristics of the 
target and the accuracy of the equipment under the expected 
experimental conditions. The absorption characteristics of the 
target indicate the thickness of the “surface” seen by an 
infrared temperature measurement device. For semitransparent 
targets, generally liquids and gases, radiation is emitted 
(absorbed) over varying depths depending on the wavelength 
of the radiation and the properties of the material. If this depth 
over which radiation is emitted is small (e.g. the total depth 
of fluid H > 1 Ik, where k is the extinction coefficient), the 
measurement of the emitted radiation is considered to be from 
the ‘ ‘surface” . The emissivity of the target must also be known 
to correct for non-ideal radiators. Most important, the ability 


l 


of the instrument to accurately measure the surface temperature 
for a given set of experimental conditions must be checked 
via calibration with other techniques. 

Commercially available infrared imaging equipment has 
been used quite extensively in the field for a variety of 
applications (refs. 4 and 5), but little documentation exists for 
its laboratory use for fluid surface temperature measurement, 
specifically, thermocapillary flow with the STDCE configura- 
tion. Moreover, there is little documentation of the absorption 
characteristics of silicone oil, the STDCE test fluid. Only one 
such reference was found (ref. 6), but it is not applicable to 
this problem. The reference data for silicone oil in solution 
cannot be used for pure silicone oil because absorption bands 
may undergo intensity or wavelength changes when the 
component is in solution (ref. 7). No reference was found 
containing data on the emissivity of silicone oil. 

In support of the STDCE an infrared thermal imaging system 
was used to measure the surface temperature of a 
thermocapillary flow under the conditions and parameters 
which are presently representative of the flight experiment. 
The absorption characteristics of silicone oil were measured 
and the absorption coefficients were calculated prior to 
acquisition of the system; these were used to determine its 
operating wavelength (region of strongest absorption) and the 
reported “surface” thickness. In order to assess the ability 
of the thermal imaging system to perform the measurements 
required by the STDCE, thermocouple measurements made 
on identical flows were compared with the IR temperature data 
and a calculated surface temperature distribution. The effective 
emissivity of silicone oil for these experimental conditions was 
also determined. The techniques and procedures used in these 
measurements are detailed, and the results and problems are 
discussed. The work described herein supported the decision 
to use this type of measurement equipment for the STDCE. 

A secondary purpose of this work was to provide guidelines 
for future users of commercially available imagers. A basic 
explanation of the types of available infrared measurement 
equipment and the performance specifications of the selected 
infrared imager are given to provide the reader with some 
background information on these types of equipment and to 
help in the selection of the device best suited for a given 
application. 

Symbols 

D diameter, cm 

d sample thickness, mm 

H depth, mm 

/ radiation intensity 

k extinction coefficient, mm -1 

Ma Marangoni number 


P pressure, Pa 

Pr Prandtl number 

p power, W 

R test cell radius, cm 

Ra Rayleigh number 

r radial position, cm 

r/R nondimensionalized radial position 
T temperature at any point on surface, °C 
z depth of thermocouple, in. 

a absorption coefficient, mm -1 

6 emissivity 

X wavelength, pm 

t transmission 

Subscripts: 
a average 

b background 

c sidewall 

r reference 

t target 

w window 

X at wavelength X 

0 initial 

1 beginning of interval 

2 end of interval 

8-12 8 to 12 pm 

10.6 10.6 pm 

Surface Tension Driven Convection 
Experiment 

Basic Description 

The basic experimental configurations for the STDCE are 
illustrated in figures 1 and 2. The figures show the two 
experimental configurations that will be studied: a constant 
temperature boundary condition (called CT experiment 
hereinafter), which will be imposed by maintaining a 
temperature difference along the free surface between the 
inside and outside walls, and a constant flux boundary 
condition (called CF experiment hereinafter), which will be 
created by heating the free surface with a radiant heater, 
thereby imposing different thermal signatures along the free 
surface. Two configurations will be studied for several 
reasons. In the CT experiment the total temperature difference 
along the free surface, which is related to the magnitude of 
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(a) Flat free surface. 

(b) Curved free surface. 

Figure 1.— Constant temperature experiment. 


the driving force (ref. 2), can be imposed and directly 
measured. In the CF experiment various heat flux distributions 
can be imposed along the free surface, and their effect on the 
flow studied (ref. 8). 

Both experimental configurations will use a cylindrical con- 
tainer 10 cm in diameter and 5 cm in height. The STDCE will 
be conducted as two discrete experiments. For the CT experi- 
ment a 1.1-cm-diam by 5-cm-high resistance heater will be 
placed in the center of the test cell. For the CF experiment 
the resistance heater will be removed, and the free surface 
heated with a radiant heater. For both the CT and the CF 
experiments the test container will be filled to different levels 
in order to create both flat and curved free surfaces. Since 
the thermocapillary driving force acts tangentially to the local 
free surface, variations in free-surface shape permit different 
driving forces to be used as an experimental parameter. A flat 
free surface greatly simplifies both the theoretical and experi- 
mental analyses and, therefore, will be used for most of the 
experiments. Although a curved free surface allows investi- 
gation of a different driving force, it is unique to a low-gravity 
environment and cannot be studied terrestrially. Consequently, 
surface temperature measurements on a curved free surface 
are not included in the present work. Test matrices, which 
include AT, Marangoni number (ref. 2), heater power, free- 
surface shape, and duration for both CT and CF experiments, 
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(a) Flat free surface, 
(b) Curved free surface. 


Figure 2.— Constant flux experiment. 


TABLE I. — CONSTANT TEMPERATURE EXPERIMENT TEST MATRIX 


Test 

number 

Temperature 
difference, 
AT, °C 

Marangoni 

number, 5 

Ma 

Heater 

power, 3 

W 

Free surface 
shape 

Duration, 

min 

8 

10 

3.9 xlO 4 

1.5 

flat 

60 

9 

25 

1.2xl0 5 

6.8 

flat 

10 

10 

60 

4.1X10 5 

32.8 

flat 

10 

11 

10 

3.2X10 4 

1.5 

curved 

10 

12 

60 

3.3 x10 s 

32.8 

curved 

10 


a Values for flat surface tests are computed numerically, and those for the curved surface are estimates. 
5 Based on average liquid depth. 


are presented in tables I and II. A more detailed explanation 
of other experimental considerations and parameters is given 
in references 1, 2, 3, and 8. 

The experimental hardware for the STDCE is presently in 
a design phase, with the basic concepts in place. The test 
chamber will be made of copper to insure good thermal 
boundary conditions at the side wall. A 10-cSt silicone oil will 
be used as the test fluid because it is insensitive to surface 
contamination (a common problem with surface tension 
experiments), it is safe, and its viscosity represents the best 
compromise between flow speed and susceptibility to free 
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TABLE n. -CONSTANT FLUX EXPERIMENT TEST MATRIX 


Test 

number 

Heating 

zone 

diameter, 

mm 

Total 

power, 

W 

Temperature 
difference, 3 
AT, °C 

Marangoni 

number, a,b 

Ma 

Free surface 
shape 

Duration, 

min 

1 

10 

0.5 

10.4 

4.2X10 4 

flat 

60 

2 

10 

3.0 

36.8 

1.5X10 5 



10 

3 

30 

3.0 

12.4 

5.0X10 4 





4 

5 

.2 

10.2 

4.1X10 4 





5 

5 

3.0 

64.7 

4.4 X10 5 





6 

5 

3.0 

64.7 

2.2X10 5 

curved 



7 

30 

3.0 

12.4 

2.5 X10 4 

curved 




“Values for the flat surface tests are computed by numerical analysis, and those for the curved surface are estimates. 
^Based on average liquid depth. 


disturbances. A planar cross section will be illuminated by 
using a 1-mm-wide helium-neon laser light sheet that is 
reflected by seed particles in the fluid, thereby allowing 
observation of the flows. Since this geometry produces an 
axisymmetric flow under most conditions, illumination of only 
one cross section will be needed. The flow data will be 
observed with an intensified video camera (which views the 
light sheet through a viewport in the copper chamber) and 
recorded with a videocassette recorder for post-flight analysis. 

Bulk flow temperatures will be measured by placing 
temperature probes in the fluid. The heating systems used will 
be a resistance heater (CT experiment) and a C0 2 laser (CF 
experiment). The infrared imager data will also be recorded 
on videotape and reduced after flight in a manner similar to 
that outlined in this paper. 

Science Requirements 

The intent of this work was to develop and calibrate an 
infrared surface temperature measurement system that will 
meet the science requirements for the STDCE. The complete 
science requirements are detailed in reference 8. The science 
requirements for the surface temperature measurement system 
are given here as a reference for the following sections. 

“To obtain information on the temperature distribution along 
the free surface a thermography technique is employed. An 
IR scanner is placed above the free surface to construct a 
thermogram of the surface, from which the temperature 
distribution is determined. The wavelength range should be 
chosen so that the measured radiation comes from the region 
as close to the surface as possible. The surface temperature 
measurement error should be within ±5% of AT and the 
spatial resolution should be less than 1 mm. To satisfy those 
requirements the technique should be carefully calibrated 
against measurement by fine thermocouples as well as against 
the numerical analysis. The response time should be less than 
0.1 seconds.” (ref. 8) 


Measurement of Absorption Coefficients 
of Silicone Oil 

Absorption Coefficients 

To understand how silicone oil emits and absorbs electro- 
magnetic radiation, knowledge of its absorption coefficients 
is necessary. Lambert’s law states that the absorbance of a 
material is directly proportional to the thickness of the material 
in the optical path. The absorption coefficient is the 
proportionality constant, which is a material property; it 
indicates how strongly the material will absorb and emit 
electromagnetic radiation. Consider a slab of an isotropic fluid 
of thickness dx with incident radiation /ox at wavelength X, 
and /ox + dl 0x the radiation leaving the slab (fig. 3). The 
decrease in radiation intensity through the slab is 

dl 0 x = ~ hloxdx ( 1 ) 

where k is the extinction coefficient, which is equal to the 
absorption coefficient a in the absence of scattering; that is, 

k(K,T,P) = a(k,T,P) (2) 

These two coefficients are, in general, functions of wavelength 
X, temperature T, and pressure P. The samples used in the 
measurement of the extinction coefficient are at the same 
pressure and temperature; therefore, 

*x = «x (3) 

Integrating equation (1) over the sample thickness d gives 

W ?d 

(dI 0X /I 0X )d\ = - k x (d*)dx (4) 

J/ x (0) Jo 



Figure 3.— Normal incident radiation on absorbing medium of thickness dx. 
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where d* is a dummy variable. Under the assumption that k 
does not vary with jc, equation (4) reduces to 

J x (</)// x (0) = exp(-M) (5) 

where I x (d)/I x ( 0) (or r) is the fraction of energy at wave- 
length X transmitted through a thickness of material d\ and 
k x is the extinction coefficient at wavelength X (ref. 9). This 
relationship is, in most cases, determined by experiment. 
Percent transmission, r X 100, was measured as a function of 
wavelength for four fluid thicknesses of 100 percent silicone 
oil. When percent transmission was plotted with respect to fluid 
thickness, the resulting curve had the form of equation (5). 
It is more convenient to take the natural log of equation (5) 
and plot -In r as a function of sample thickness d. This yields 
a straight line whose slope is equal to the absorption coeffi- 
cient. Since infrared temperature devices generally operate in 
a specific wavelength range because of atmospheric trans- 
mission losses in certain regions, an average transmission T a 
was calculated for a wavelength range X! to X 2 , 


>X 2 

rd\ 



The integral was evaluated numerically by using Simpson’s 
rule (ref. 10). The wavelength range for the calculated average 
transmission was determined by the operational wavelength 
ranges of available IR equipment and the regions where 
silicone oil best absorbs electromagnetic radiation. Since the 
C0 2 laser radiation determines the thermocapillary flow 
driving force, the absorption coefficient is also calculated at 
the C0 2 laser wavelength (approximately 10.6 /an). 

Apparatus and Procedure 

A Perkin-Elmer Model 780 spectrophotometer was used to 
measure the percent transmission 7 as a function of wave 
number, the reciprocal of wavelength. The wave number range 
covered by the spectrophotometer is 4000 to 500 cm -1 
(2.5 to 20 pm). The spectrophotometer has two optical paths— a 
measurement path and a reference path. A sample holder 
(fig. 4) consisting of two sodium chloride (NaCl) windows 
(1.5 in. high, 0.75 in. wide, and 0.25 in. thick) separated by 
a spacer of variable thickness is placed in the measurement 
path, and a dummy sample holder, in the reference path. 
Sodium chloride is used for the window material because it 
has a flat transmission curve in the infrared. The sample of 
fluid is injected into the cavity between the two windows. The 
dummy sample of a single 0.5-in. -thick NaCl window is placed 
in the reference beam to account for the transmission loss 
through the NaCl windows in the sample holder. 



Figure 4.— Fluid sample holder for Perkin-Elmer Model 780 spectrophometer. 

All control commands for the spectrophotometer are executed 
through a personal computer (PC) interfaced with the spectro- 
photometer. During operation the spectrophotometer emits a 
single beam which is split into two beams of constant intensity 
but varying wave number. After the beams pass through the 
respective samples, the intensity of each beam is measured. 
The two measurements are compared, and the result, percent 
7, is calculated and plotted as a function of wave number. The 
data and associated graph are then stored on a floppy disk for 
later analysis. A block diagram of the apparatus is shown in 
figure 5. 

Four tests were conducted at 25 °C, with different sample 
thicknesses: 0.025, 0.050, 0.100, and 0.200 mm. After 
smoothing the data to eliminate small peaks caused by noise 
in the spectrophotometer, the percent 7 as a function of wave 
number was extracted via a user program on the PC. 

Absorption coefficient in the 8 - to 12-pm region .— To 
calculate the percent r a , a user program was written to extract 
data at intervals of 10 wave numbers in the region 1250 cm -1 
to 830 cm -1 (8 to 12 /mi). The program, using Simpson’s 
rule, calculates the area under the curve of percent 7 as a 
function of wave number; then it divides this result by the 
interval length (420 cm -1 ) to give percent r a . This procedure 
was repeated for each of the four sample thicknesses. The 
average percent transmission was plotted with respect to 
sample thickness; this curve is the empirical relationship 
described by equation (5). These data were plotted again as 
—In 7 with respect to sample thickness, and a linear least 
squares regression (ref. 11) yielded a, the slope of the line. 
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Figure 5.— Spectrophotometer block diagram. 


Absorption coefficient at 10,6 pm.— The expected output 
wavelength of the C0 2 laser used by the STDCE is 10.6 pm 
(943 cm" 1 ). The values of percent r at wave numbers 
940 cm -1 (10.64 pm) and 950 cm -1 (10.52 pm) were 
averaged for each of the four sample thicknesses, and from 
the procedure described previously, the absorption coefficient 
at this wavelength was determined. 

Results and Discussion of Absorption 
Coefficient Measurements 

The absorption coefficient results are summarized in tables 
m and IV. The raw data of percent rasa function of wave 
number for the four sample thicknesses are presented in table V. 
Figure 6 is a plot of percent r as a function of wavelength 
for 100 percent silicone oil in the 4000- to 500-cm _1 (2.5- 
to 20- pm) region with fluid thicknesses of 0. 100, 0.025, and 
0.200 mm. The data for the sample thickness of 0.050 mm 
were not stored on a disk and, therefore, are not available in 


TABLE in.-AVERAGE TRANSMISSION r a 
AND -In T a AS FUNCTIONS OF 
SAMPLE THICKNESS 


Sample 
thickness, 
d, mm 

Wavelength range, X 

/xm 

8 to 12 

10.6 


-In r a 

7 a 

-In 7 a 

0.000 

1.000 

0.00 

1.000 

0.00 

.025 

.225 

1.41 

.657 

.42 

.050 

.130 

2.05 

.426 

.85 

.100 

.045 

3.11 

.178 

1.73 

.200 

.003 

5.74 

.027 

3.60 


TABLE IV. — ABSORPTION COEFFICIENTS 
FROM LEAST SQUARES REGRESSION 


Wavelength range, 

Absorption coefficient. 

X, /xm 

a, mm -1 

8 to 12 

27.0 

10.6 

18.0 


TABLE V. — TRANSMISSION AS A FUNCTION 
OF WAVE NUMBER FOR FOUR 
SAMPLE THICKNESSES 
(RAW DATA) 


Wave 

number, 

1/X, 

cm' 1 

Sample thickness, 
d, mm 

0.025 

0.050 

0.100 

0.200 

1250 

14.25 

4.18 

0.25 

-0.54 

1240 

38.81 

17.95 

4.70 

-.09 

1230 

59.75 

34.69 

12.40 

1.59 

1220 

67.22 

45.05 

19.68 

3.62 

1210 

68.24 

47.06 

21.85 

4.15 

1200 

65.49 

42.78 

18.57 

2.80 

1190 

58.49 

33.68 

11.81 

0.99 

1180 

47.69 

22.75 

5.65 

-.05 

1170 

35.02 

12.73 

1.88 

-.31 

1160 

22.92 

5.57 

.33 

-.29 

1150 

12.31 

1.42 

-.31 

-.38 

1140 

4.61 

-.14 

-.63 

-.55 

1130 

.51 

-.49 

-.65 

-.41 

1120 

-.39 

-.46 

-.54 

-.19 

1110 

-.43 

-.44 

-.42 

-.36 

1100 

-.48 

-.48 

-.40 

-.57 

1090 

-.77 

— .46 

-.35 

-.68 

1080 

-.76 

-.42 

-.02 

— .46 

1070 

-.62 

.01 

.14 

-.19 

1060 

-.12 

.46 

.12 

-.05 

1050 

.07 

.37 

-.09 

-.03 

1040 

.04 

-.11 

-.09 

.02 

1030 

-.27 

-.20 

-.06 

-.06 

1020 

-.77 

-.30 ! 

-.05 

-.09 

1010 

-.34 

-.47 

-.17 

-.28 

1000 

5.39 

.06 

-.22 

-.33 

990 

18.03 

3.98 

.13 

-.42 

980 

33.71 

11.72 ! 

1.60 

-.17 

970 

47.21 

22.13 

5.02 

.09 

960 

57.24 

32.39 

10.54 

.91 

950 

64.10 

40.51 

16.12 

2.26 

940 

67.19 

44.61 

19.49 

3.22 

930 

66.19 

43.26 

18.59 

2.75 

920 

61.40 

37.20 

13.93 

1.35 

910 

52.12 

26.71 

7.56 

.17 

900 

37.00 

14.31 

2.29 

.56 

890 

19.79 

4.40 

-.33 

-.74 

880 

7.70 

-.0 

-.78 

-.62 

870 

3.18 

-.66 

-.29 

-.56 

860 

2.00 

-.40 

.11 

-.64 

850 

.48 

-.42 

.15 

-.60 

840 

-.39 

-.42 

-.28 

-.46 

830 

-.27 

-.40 

-.56 

-.52 
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TRANSMISSION, T, PERCENT 



WAVELENGTH, A, MM 

(a) Thickness = 0.025 mm. 

(b) Thickness = 0.100 mm. 

(c) Thickness = 0.200 mm. 

Figure 6.— Transmission as a function of wave number for three silicone oil sample thicknesses. 
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graphic form. These plots show that the region of strongest 
absorption was between 1250 and 700 cm -1 (8 and 14 pm). 
This result combined with the availability of an imager which 
operates in the 8- to 14-/rm region determined the region to 
be averaged. The percent r a was calculated for the region 
between 1250 and 830 cm" 1 (8 and 12 pm) because the 
commercial imagers that were available at the time operated 
in this region. Adding the contribution of the 12- to \A-pm 
region would not have a profound impact on the “surface” 
as seen by the imager since the imager detector response 
decreases sharply in this region (as will be discussed later), 
thus de-emphasizing its contribution. 

Table III lists the T a of sample thicknesses for both wave- 
length ranges; these data are shown graphically in figure 7. 
The results of the regression on both data sets are listed in 
table IV and are plotted in figure 8. In figure 7 the reduction 
in transmission of radiation by a factor of 1/e (or 63.2 percent 
absorbed) in the 8- to 12-pm region occurs at a sample thick- 
ness of approximately 0.012 mm. Similarly, for radiation at 
10.6 pm the sample thickness is 0.060 mm. This thickness 
is commonly known as the mean penetration length and is 
larger for the data at 10.6 /mi because there is a minor trans- 
mission window around that wavelength. The transmission in 
the 8- to 12-pm region at a sample thickness of 0.200 mm 
is nearly zero. Therefore the “surface” as seen by an imager 
is approximately 0.200 mm thick. A comparison of the results 
shown graphically in figure 6 and tabularly in table HI indicates 
that the C0 2 laser radiation penetrates deeper than the surface 
the imager sees. This is advantageous since, if the situation 
were opposite, the temperature indicated by an infrared 
temperature measurement device would be an average 
temperature of the region where a temperature gradient exists 
due to the C0 2 laser heating. 


i.o? 

i 
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MM 

8 TO 12 
10.6 


TRANSMISSION = 1/e 
(63 PERCENT ABSORBED) 


0 .02 .04 .06 .08 .10 .12 .14 .16 .18 .20 

SAMPLE THICKNESS, d, mm 


Figure 7.— Average transmission of silicone oil as a function of sample 
thickness d. 
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Figure 8.— Plot of -In of average transmission as a function of sample 
thickness d. Slope of line is absorption coefficient a. 


Surface Temperature Measurements 

Apparatus 

Infrared temperature measurement equipment. —Presently 
there are several types of noncontact infrared temperature 
measurement devices of differing levels of sophistication: spot 
radiometers, qualitative infrared imagers, and quantitative 
infrared imagers (ref. 12). The basic operation of the three 
types is similar. All targets radiate energy whose intensity is 
a function of temperature and wavelength. For a black-body 
the functional relationship is the Planck function, and for a 
gray-body the functional relationship is proportional to the 
Planck function, where the proportionality constant is called 
the emissivity. The emitted energy is received by a detector 
in the device that produces an output signal related to the 
emitted intensity. 

Spot radiometers (also called infrared pyrometers), in 
general, use pyroelectric detectors, which develop an electrical 
surface charge when the surface temperature is changed. These 


instruments are small and easily portable but have two major 
drawbacks— a single point measurement and a large response 
time (seconds). 

Both qualitative and quantitative infrared imagers (also called 
thermographic devices) use quantum detectors (photodetectors), 
such as mercury cadmium telluride semiconductors in which 
the electrons are raised from nonconducting to conducting 
states by the energy contained in photons arriving from the 
target. The result is an electrical charge proportional to the 
intensity of the incoming radiation. Imagers produce a full two- 
dimensional temperature field by scanning a detector field of 
view, called an instantaneous field of view (IFOV), either over 
the entire field of view of the imager (scanning imagers) or 
by staring at the target with an array of detectors, called a focal 
plane array. The output is either displayed at 1/30 sec (for 
TV-compatible systems) on a video monitor or digitized and 
stored on a disk. Imagers overcome the previously mentioned 
problems of spot radiometers, but they introduce other 
constraints. Quantum detectors require cooling (generally to 
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liquid nitrogen temperatures) to achieve sensitivity to infrared 
wavelengths, and they are sophisticated, delicate instruments, 
which reduces their portability. 

The distinction between quantitative and qualitative imagers 
is determined by whether the instrument has the ability to 
measure absolute temperatures or merely display intensities 
(e.g. night vision equipment). The quantitative imager has the 
ability to extract accurate temperature data from the image, 
whereas the qualitative imager provides the ability to detect 
and recognize active radiation emitted from a target (as 
opposed to visible radiation, which is reflected). 

Selection criteria .— The class of infrared imagers that is 
most suited for the STDCE, based on the aforementioned 
characteristics, is a quantitative infrared imager. As stated 
previously, the STDCE requires a full field noncontact tempera- 
ture measurement device that has a response time of less than 
0. 1 sec, a spatial resolution of 1 mm, and the ability to extract 
qualitative surface temperature data from the image. The 
commercially available infrared imagers (henceforth, imager 
refers to a scanning imager) are typically offered with spectral 
sensitivities in three wavelength ranges: short band (1 to 
5 fim), long band (8 to 14 fim ), and broad band (3 to 14 pm). 
The wavelength range of the imager is decided by two factors— 
the absorption characteristics of the target and the expected 
temperatures of the target. The expected STDCE surface 
temperatures are less than 100 °C above normal ambient 
temperatures (25 °C). Since the peak wavelength of a black- 
body at these temperatures is approximately 10 pm, an imager 
with a peak detector sensitivity near this point is desirable to 
achieve the highest signal-to-noise ratio. The area of strongest 
absorption for silicone oil is between 8 and 14 pm; therefore 
a long band imager is best suited for use with the STDCE. 

A long band infrared imager and image processing system 
were purchased from Inframetrics, Inc. The state-of-the-art 
measurement performance of the Inframetrics Model 600 
infrared imager and its user features make it an excellent choice 
for this particular application. 

Inframetrics thermal image processing system — The target 
absorption coefficient must be known and the performance of 
the imaging system understood before the measurements can 
be interpreted. Specifications such as detector type, band pass, 
thermal sensitivity, information rate, image resolution, 
dynamic range, and calculation of absolute temperatures are 
the factors that determine how the imager measures radiation 
from, and calculates the temperature of, a surface. 

The basic operational concept of the Inframetrics Model 600 
and image processing system is presented herein to provide 
a general understanding of the equipment. Detailed descrip- 
tions of each of the performance specifications are presented 
in the appendix for readers interested in understanding the 
complete operation of the Inframetrics system. 

Model 600 thermal infrared imager : The Inframetrics 
infrared imager (henceforth called the Model 600) is a 
scanning-type imager. An instantaneous field of view (IFOV) 
is scanned over the field of view (FOV) of the scanner. 


Radiation from the IFOV is focused on the mercury cadmium 
telluride detector as the IFOV is scanned over the FOV at 60 
frames per sec, thereby constructing the thermographic image. 
The image is digitized into discrete gray levels, and the 
intensity of each pixel is converted in real time to a temperature 
via factory supplied calibration curves. The image is converted 
to a standard TV signal (EIA RS-170 or NTSC composite) 
and displayed on a standard TV monitor in false color 
proportional to temperature. User inputs to the system are 
target emissivity (effective emissivity) and background 
temperature; these account for nonblack-body radiation and 
background reflection radiation from the target. The spatial 
resolution of the Model 600 is 0.8 mm (for a working distance 
from the scanner to the target of 40 cm), and the response 
time is 0.033 sec. Both satisfy the science requirements. 

The Model 600 has many features available to the user. The 
temperature ranges available are 5, 10, 20, 50, 100, and 
200 °C (normal range) and 500 and 1000 °C (extended range). 
These temperature ranges are called the temperature window. 
The center temperature of the window is user selectable 
thereby allowing these two variables to be combined to 
maintain full sensitivity— that is, the smallest number of, or 
fraction of, degrees per gray level while keeping the 
temperature of all desired objects within the minimum and 
maximum selected temperature. 

As the image is displayed on the TV monitor the user can 
select any one of three measurement modes: point, area, and 
line scan. In point mode the temperature of a single picture 
element (pixel) is reported when a set of cross hairs is 
positioned on that pixel. Line scan mode allows the user to 
select a cross section of an image and display the temperature 
profile of that cross section. These modes, as well as target 
emissivity, background temperature, and a variety of other 
housekeeping functions, are selected from the keypad on the 
control unit. 

Thermoteknix thermal image processing system : The combin- 
ation of the Model 600 scanner, electronics and control unit, 
and monitor is a stand-alone, quantitative thermal imaging 
system capable of solving the radiometric equation in real time. 
The system can be made even more powerful with the addition 
of the Thermoteknix thermal image processing system. The 
Thermoteknix system is an IBM PC/AT-based dedicated image 
processing and measurement system that includes all features 
available on the standard Model 600 (point, line, and area 
modes) plus a plethora of other functions. These functions 
include continuous 8-bit digitization, digital disk image storage, 
three point measurement modes, thermal gradient, image 
magnification, thermal magnification, image averaging, and 
color hardcopy. Since the image processing system contains 
all the necessary calibration curves to convert intensity levels 
from the image into absolute temperatures, the Model 600 
electronics and control unit is used only for construction of 
the image, scan conversion to television format, temperature 
range and center temperature selection, and some house- 
keeping details. 
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In addition to the Model 600 and Thermoteknix image 
processing system (shown in fig. 9), the experimental hardware 
consists of the test cell, test fluid, heating systems, and 
thermocouple instrumentation. All hardware is assembled on 
a Melles Griot breadboard (fig. 10). 

Test cell and fluid.— The test cell used for both the CT and 
CF experiments was a 10.2-cm-diam by 6.2-cm-high copper 
cylinder with a 0. 16-cm wall thickness, which is the approxi- 
mate size of the test cell required in the STDCE. A 0. 16-cm 
copper plate was soldered to the end of the copper cylinder 
to provide a bottom. In order to maintain a constant tempera- 
ture boundary condition at the side and bottom walls, constant 
temperature water was circulated through a 1.0-cm-diam 
copper tube soldered to the outside diameter and bottom of 


the test cell (fig. 1 1). A Hotpack Corp. Model 603340 constant 
temperature recirculating bath circulated constant temperature 
water through the tubing at approximately 12 liters/min with 
a temperature control of ±0.01 °C. A plexiglas shield was 
placed over the test cell to reduce reflected background 
radiation. The fluid used for these experiments and with the 
STDCE was Dow Corning 200 fluid (silicone oil) of 10-cSt 
viscosity. 

Heating systems— The external heating system for the CF 
experiment was an 8-W RF-excited C0 2 laser manufactured 
by Laakman Electro-optics (now Melles Griot Laser Products). 
The laser delivered a constant energy flux to the fluid surface 
at normal incidence. Since the C0 2 laser beam is outside 
visible wavelength region, a 5-mW Spectra Physics Model 159 
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helium-neon laser facilitated alignment of the beam as shown 
in figure 12. The He-Ne laser was combined with the C0 2 
beam by means of a beam combiner manufactured by Two- 
Six, Inc. This optic reflected the He-Ne beam while trans- 
mitting the C0 2 beam, thereby allowing concentric alignment 
of the two beams. Since the two beams followed identical 
paths, the optics were aligned only when the C0 2 laser was 
off and the He-Ne laser on, thereby reducing the possibility 
of injury caused by the more dangerous C0 2 laser. The laser 
power was read in watts directly from a Coherent, Inc. , Model 
201 power meter. Because the C0 2 laser required such a 
large power input (225 W) to deliver a small output, there 
was a considerable amount of waste heat that was removed 
via a water-cooled baseplate. The output stability of the laser 
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Figure 12.— Optical paths. 


was dependent on the thermal control of the laser cavity. Tap 
water provided a coolant that was stable enough to maintain 
the laser output to approximately a 5-percent variation. 
Because of the beam geometry, the C0 2 laser could produce 
a 5-mm-diam beam without any lenses in its optical path. 

A cartridge heater press-fit into a 1 . 1-cm-o.d. by 5-cm-high 
copper shell heated the fluid for the CT experiment. The 
cartridge is essentially a resistance that produces heat when 
an electrical current is passed through the circuit. For the CT 
experiment the center of the test cell was located; then the 
cartridge heater inside the copper shell was fastened to the 
bottom wall with epoxy cement (fig. 11). This configuration 
yielded a constant temperature boundary condition between 
the heater and the side wall. A variable AC power supply 
provided the power source for the heater. The voltage and 
current were measured with a John Fluke Mfg. Co., Inc., 
Model 8840 A multimeter. 

Thermocouple instrumentation . —Thermocouples measured 
the flow, wall, and background temperatures. The background 
temperature was measured with a type-E thermocouple 
mounted on the plexiglas radiation shield above the test 
chamber. This location was chosen because the background 
radiation that was reflected from the fluid surface and seen 
by the IR imager originated from this surface. Therefore, this 
temperature was the most accurate measure of the background 
radiation. During the CT experiment, type-T thermocouples 
were mounted on the heater and side walls, whereas during the 
CF experiment only the sidewall thermocouple was used 
because the cartridge heater was removed. For the thermocouple 
surface temperature measurement, a type-E thermocouple with 
a 0.001 -in. bead was attached to a J-shaped probe positioned 
by using a three-axis linear translation stage so that the probe 
did not penetrate the surface in the vicinity of the bead 
(fig. 13). The thermocouples’ emf was read on the Fluke 
multimeter with electronic cold junction compensation per- 
formed by Omega Engineering Corp.-CJ (type-E and type-T) 
cold junction compensators. 



Figure 13.— J-shaped thermocouple probe and measurement location. 
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Procedure 

Methodology .— With the knowledge of the absorption 
characteristics of silicone oil and an understanding of the 
operating and performance characteristics of the Model 600, 
surface temperature measurements were made and interpreted. 
The understanding of these characteristics provided the 
background to explain the sources of error in the infrared 
surface temperature measurements. 

Two other methods were used as standards with which to 
compare the data from the Model 600. For the CT experiment 
(internal heating), thermocouple measurements and a 
calculated surface temperature distribution were used. For the 
CF experiment (external heating), only the calculated surface 
temperatures were used, for reasons which will be explained 
in the section Thermocouple Measurements and Calculated 
Temperature Distribution. 

Constant temperature experiment.— The test cell was 
washed with soap and hot water and leveled on the breadboard. 
A barrier coating (3M FC-723) was applied to the heater edge 
to ensure that the fluid would not spill over the edge (fig. 11). 
The test cell was then filled to the edge of the heater wall (5 cm 
from the bottom wall) with silicone oil. Before data were taken, 
the distance between the edge of the cartridge heater and side 
wall at the azimuthal angle 6 (fig. 14), where the measurements 
were to be taken, was measured with a dial caliper. This length 
plus half the heater diameter was used as the test cell radius 
R to nondimensionalize the radial position (i.e., r/R ). The 
three-dimensional positioner was adjusted such that one axis 
was parallel to the desired radius. The thermocouple was 
traversed along and normal to this radius for each measurement 
by adjusting the translation stage micrometers, which are 
calibrated to 0.01 mm. The pixel locations in the IR image 
that corresponded to the center of the heater and the side wall 
were determined. The distance between these locations was 


z 



used to nondimensionalize the radial position for the IR Imger 
data. The intermediate pixels were addressed, and the 
temperature distribution as a function of pixel location was 
extracted. 

Power applied to the cartridge heater and constant 
temperature water (23 °C setpoint) circulating around the test 
cell initiated the experiment. The experiment was run for about 
2 hr so that the surface temperature could reach steady state 
for comparison with the calculated surface temperature 
distribution, which was computed assuming a steady-state 
model. 

At the start of the run the IR image was stored. The image 
used for extraction of temperatures was the average of 10 
frames acquired at 1 sec intervals. This reduced the noise in 
the image by a factor of the square root of the number of 
frames averaged. Immediately following the storage of the 
image, the surface temperature was measured by traversing 
the thermocouple along the measurement radius at intervals 
of 0.5-mm from r/R = 0.1 to 0.24 (near the heater); 1 mm 
from r/R = 0.24 to 0.34; 5 mm from r/R = 0.34 to 0.94; and 
1 mm from r/R = 0.94 to 1.00. The grid size was changed 
to accommodate a sharp temperature gradient near the walls. 
The position of the thermocouple in relation to the surface is 
extremely critical. The accepted method for positioning the 
thermocouple is to traverse the probe, starting from the bulk, 
normal to the surface until the bead begins to distort the 
surface. At this position the reading was taken (fig. 13). The 
position was readjusted for each measurement point. Following 
the acquisition of surface temperature data, the side wall, 
heater, and background temperatures were again measured to 
determine whether the experimental conditions had changed. 

Constant flux experiment.— Between the CT and CF 
experiments, the test cell was disassembled (the cartridge 
heater removed), washed with soap and hot water, releveled, 
and secured to the breadboard. The test cell was filled with 
clean silicone oil to a level of 5.5 cm from the bottom. The 
procedure for aligning the thermocouple, outlined previously, 
was repeated with the following changes. Since there was no 
inner wall (cartridge heater), the center of the test cell was 
determined by the intersection of two diameters at right angles. 
Again, the distance from this point to the side wall, along the 
radius where the measurement would be made, was measured. 
The pixel location of the side wall, along the measurement 
radius, in an IR image of the test cell was determined. 

The C0 2 laser beam was aimed at the center of the test cell 
(as determined above) at nearly normal incidence. Fine 
adjustments to the beam position were made by rotating the 
turning mirror that directed the beam toward the fluid surface. 
The beam was adjusted until the output of the thermocouple 
placed at the center of the test cell was at a maximum, thereby 
indicating that the maximum intensity (center) of the beam 
coincided with the center of the test cell. The pixel location 
reporting the highest temperature on the IR image was 
recorded as r = 0. 
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After powering up the laser, the output power was stabilized 
(0.5 hr) with the test cell cooling water at the 23 °C setpoint. 
To initiate the experiment the beam block was removed, 
thereby allowing the laser to heat the surface. Again the 
experiment was run for about 2 hr to allow a steady-state 
condition to be reached. The background and wall temper- 
atures were recorded at the beginning and end of the data 
acquisition. The power meter was momentarily placed in the 
beam path to record the laser power. After removing the power 
meter and waiting 5 min, when the thermocouple output was 
again steady, the IR image was averaged and recorded in the 
same manner as before. The thermocouple was traversed 
across the surface, and measurements were taken at radial 
locations identical to the CT experiment, except the measure- 
ments extended to the center at intervals of 0.5 mm. 

Data reduction .— The thermocouple data reduction 
consisted of converting the measured voltage to temperature 
with a basic computer program that incorporates the Omega 
Engineering Corp. temperature-to- voltage chart (type E). The 
program linearly interpolated between values in the table to 
determine the temperature of the corresponding measured 
voltage. The results where plotted as a temperature difference 
(T — T c ), where T c is the thermocouple-measured side wall 
temperature, as a function of radial position rfR. 

The surface temperature data from the IR imager was 
extracted from the image stored on the hard disk of the image 
processing computer. Since the end points of the measurement 
radius were determined, a nondimensional radial position, in 
terms of screen pixels, was determined. After the background 
temperature and target emissivity were input, the surface 
temperatures at specified pixel locations were extracted from 
the image by the point measurement mode. The pixel locations 
sampled were determined by using a grid size similar to that 
used with the thermocouple measurements. A variable grid 
size was used because of the sharp gradient near rlR = 0 and 
r/R = 1 . These data were plotted as a temperature difference 
(T - T c ) where T c was the side wall temperature measured 
by a thermocouple. This difference was plotted with respect 
to radial position r/R on the same graph with the calculated 
surface temperature distributions. 

Results and Discussion 

Effective Emissivity of Silicone Oil 

The user inputs to the Model 600 are the background temper- 
ature and the target emissivity. Arriving at the background 
temperature measurement was not a problem because it can 
easily be measured with a variety of devices. Determining the 
target emissivity, on the other hand, was not as straight- 
forward. Since the emissivity of silicone oil was not to be found 
in the literature, the effective emissivity was determined 
empirically. This value is called the effective value because 
it was not measured but implied by comparison of the IR 


imager surface temperature data with the thermocouple and 
calculated surface temperature distributions. This implied 
emissivity, or effective emissivity, is the effective average^ 
spectral directional emissivity. The spectral characteristic was 
determined by the operating wavelength of the imager (8 to 
14 /on), whereas the directional characteristic was determined 
by the angle at which the imager viewed the surface. In order 
to accommodate the optical path of the C0 2 laser (fig. 12), 
the imager must view the target at a small angle from the unit 
normal. In reference 9 the directional spectral emissivity is 
predicted from electromagnetic theory as a function of emitted 
ray angle /? measured from the unit normal. However, these 
results indicate that the directional spectral emissivity does not 
vary until 0 becomes larger than 45° (ref. 9). Therefore, 
viewing the surface at an angle of approximately 10° should 
not affect the reported surface temperature. 

Thermocouple Measurements and Calculated Surface 
Temperature Distribution 

Before thermocouple measurements were made for the CT 
and CF experiments, Kamotani (ref. 3) calculated the 
temperature distribution using a numerical model developed 
to supplement the STDCE and to assist in its design. The 
thermocouple measurements made in the CT experiment by 
using the aforementioned procedures agree very well with the 
calculated surface temperature distribution. 

A major problem was encountered with the thermocouple 
measurements made in the CF experiment. Many of the 
thermocouple surface temperature measurements made were 
not repeatable. The problem was caused by attempting to 
measure the surface temperature under the laser beam by using 
the positioning technique previously described. The thermo- 
couple was heated directly by the laser, which resulted in a 
distorted measurement. This direct heating made it impossible 
to determine if the thermocouple was reporting the flow 
temperature or the temperature resulting from direct heating. 
Therefore, the temperature reported by the thermocouple 
fluctuated greatly depending on its placement in relation to 
the surface and heating depth. Figure 15 shows the surface 
temperature distribution measured by the thermocouple at 
different vertical (z) locations as a function of r/R. The 
thermocouple was traversed across the surface five times at 
different vertical locations, given in thousandths of an inch. 
These numbers are a measure of the thermocouple position 
relative to an arbitrary zero. In order to determine the relation- 
ship between the thermal boundary layer and the heating depth, 
measurements were made (fig. 16) by traversing the thermo- 
couple in the z-direction starting from the bulk fluid. The line 
between the thermal boundary layer and the direct heating 
region is not well defined, but the approximate boundary layer 
thickness is 0. 15 mm (0.006 in.) Recall that the depth at which 
the thermocouple should be placed to avoid direct heating is 
0.265 mm (i.e., the depth at which 99 percent of the heater 
radiation is absorbed). The thermocouple could not be placed 
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thermocouple locations for constant flux experiment. Zero location is 
arbitratry, /? = 2.7 W, and beam diameter - 5 mm. 
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both beyond this absorption depth and in the thermal boundary 
layer. For these reasons the thermocouple surface temperature 
measurements were not used to assess the accuracy of the IR 
imager measurements. 

Constant Temperature Experiment 

In all experimental runs the difference between the beginning 
and ending wall, background, and heater temperatures was 
less than 0.2 °C for each. As mentioned previously, the effec- 


tive emissivity was determined empirically by comparison of 
imager data with thermocouple and calculated surface temper- 
ature distributions. The effective emissivity was varied until 
the best agreement between the IR imager data and the data 
from other techniques was achieved. On the basis of experience 
with a spot radiometer and qualitative experience with the 
Model 600, the effective emissivity was placed near 0.9. 
Figure 17 is a plot of temperature difference (T — T c ) as a 
function of radial position rlR for the three techniques, with 
three values of effective emissivity for the imager technique. 
The agreement between the thermocouple measurements and 
the calculated surface temperature distribution is very close, 
which is a confirmation of the accuracy of the numerical model. 
The values of effective emissivity plotted are e = 0.95, 0.90, 
and 0.85. Of these, the value of 0.90 resulted in the best com- 
parison over the temperature range produced under the experi- 
mental conditions (input power to the cartridge heater = 17.7 
W and wall temperature = 23.0 °C). From these data it was 
evident that different values of effective emissivity produce 
a better comparison in certain regions (temperature ranges). 
Effective emissivity, then, is dependent on temperature. In 
regions of higher temperature the value e =0.95 produced the 
best comparison, and in regions of lower temperature the value 
€ = 0.85 produced the best comparison. 

The data presented in figures 18 and 19 show two different 
experimental conditions (input power levels). The choice of 
experimental conditions was decided by the temperatures 
expected to be encountered in the STDCE. Figure 18 is 
identical to figure 17 except that the only value of effective 
emissivity plotted is e = 0.90. A photograph of the actual 
image from which the IR temperature data were taken is shown 
in figure 20. The temperature (in degrees C) is proportional 
to the color, as indicated in the scale on the left side of the 
photograph. This condition represents the probable maximum 
temperature difference expected in the STDCE for the CT 
experiment. Figure 19 shows the CT experiment with a power 
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Figure 17.— Comparison of thermocouple and calculated surface temperature 
distributions with those of IR imager at three emissivities for constant 
temperature experiment with p = 17.7 W and T c - 23.0 °C. Numerical 
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Figure 19.— Comparison of the numerical, thermcouple, and IR imager surface 
temperature distributions for the constant temperature experiment with 
€ = 0.90, p = 5.0 W, and T c = 22.8 °C. The numerical calculation used 
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Figure 20.— Constant temperature experiment IR image for T = 43.4 °C, R = 50, e = 0.90, and p = 17.7 W. 



input of 5.0 W. The agreement among the three techniques 
for the 17.7- and 5.0-W power levels is ±5 percent of the 
AT between the heater and side wall in the region r/R = 0.125 
to 0.90. Outside this region (or near the side walls) the 
measurement errors increase substantially to approximately 
± 10 percent of the AT. These errors can be explained by using 
the imager slit response function (SRF). (See the appendix for 
a detailed discussion of the slit response function.) Since there 
is such a sharp temperature gradient near the side wall, the 
imager does not follow exactly the nearly step change in 
temperature. The effect of the slit response function is given by 

^apparent = SRF X T 0 bj ect + (1 — SRF) X T surround 


Since this gradient exists over a smaller distance than a 
5-IFOV length, some measurement error is expected. To 
compound the problem, the emissivity of the copper test cell 
is very low (approximately 0.10), which reduces the radiation 
intensity seen by the imager (i.e., reports a low temperature) 
and exaggerates the displayed temperature gradient. The error 
near the heater wall is confined to a more narrow region 
( r/R = 0.1 to 0.17) but is similar in magnitude to the error 
at the side wall. In this case the lower emissivity surface is 
at a higher temperature than the surroundings; the increase 
in the radiation intensity due to higher temperature is offset 
by the loss in intensity due to the low emissivity, and thus, 
the emissivity effect is less pronounced. Applying a high 
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emissivity coating to the heater and copper chamber surfaces 
might offer a partial solution to this problem. This coating 
would aid in the resolution of the temperature gradient. Since 
the height of the test chamber is 6.2 cm but it is only filled 
to 5 cm, a meniscus is formed at the side wall. There is no 
meniscus at the heater wall because the fluid level is equal 
to the level of the heater edge (fig. 21). Silicone oil exhibits 
a near-zero contact angle on copper, and therefore, the emitted 
ray angle will range from 0° to 90°. The presence of this 
meniscus at the side wall may influence the reported 
temperature because of the previously described angular 
dependence of the emissivity. 

Because a thermocouple measurement was taken at the 
heater wall, we know the inside endpoint of the surface 
temperature distribution and can, therefore, extrapolate the 
imager data to this point. The adoption of the above solutions 
should help decrease the measurement errors at the walls. 

Constant Flux Experiment 

Figure 22 shows a typical surface temperature distribution 
under the constant heat flux condition. The power level and 
beam diameter were set by the STDCE requirements. The most 


difficult CF experiment measurement case for the STDCE will 
be at the highest power level and smallest beam diameter (3 W 
and 5 mm), because these parameters produce the sharpest 
temperature gradients, which are the most difficult to resolve. 
For figure 22 the laser power was 2.7 W (3.0 W measured 
minus 10 percent reflected), and the wall temperature was 
22.7 °C. Analysis of this plot shows that the IR imager 
distribution is within =t 5 percent of the numerical calculation 
distribution with the exception of several points in the region 
r/R = 0.05 to 0.15. Again, measurement errors are present 
at the side wall and are attributed to the problems mentioned 
in the CT experiment. There is some error at r/R = 0 caused 
by the imager reporting the average temperature of a layer 
0.2 mm thick. The boundary layer temperature changes greatly 
over this layer, so the imager was expected to report a lower 
temperature than the numerical calculation because it was 
averaging the temperature of the boundary layer and the cooler 
bulk fluid. Figure 23 is an image of the thermocapillary flow 
under these same conditions, but it is not the image used for 
data reduction. Figure 24, which is the identical image of 
figure 23 shown in relief, gives a qualitative picture of the 
temperature distribution. 
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Figure 21.— Meniscus at edge of test cell. 
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Figure 22.— Comparison of numerical and IR imager surface temperature 
distribution for constant flux experiment with 6=0.90, p = 2.1 W, 
T c —221 °C, and beam diameter =5 mm. Numerical calculation used 
Ma=1.27x 10 5 and Ra=4.87xl0 7 . 
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Concluding Remarks 

The work described herein was performed as support for 
the Surface Tension Driven Convection Experiment. The 
infrared imaging system was investigated to determine the 
feasability of using this type of measurement device for the 
STDCE. The following conclusions are drawn from this work: 

1 . An infrared imager operating in the 8- to 12- pm region, 
with a spectral bandpass similar to that of the Inframetrics 
Model 600, measured the temperature of a layer of silicone 
that was approximately 0.200 mm thick. 

2. Using a calculated absorption coefficient for silicone oil 
allowed the C0 2 laser radiation to be determined as 99 
percent absorbed over a depth of 0.256 mm of silicone oil. 

3. The effective spectral directional emissivity of silicone 
oil for these experiments was 0.90. 

4. The Model 600 infrared imager measured the surface 
temperature of the thermocapillary flow for the constant 
temperature experiment to within ± 5 percent of the overall 
AT (except for regions near the heater wall and side wall), 
based on the comparisons with thermocouple measurements 
and numerical calculations. In the excepted regions, measure- 
ment errors were as high as ± 10 percent of the overall AT. 


5. For the constant flux experiment, the Model 600 
measured the surface temperature of the capillary flow to 
within ±5 percent of the overall AT except for the region near 
the side wall and at r/R = 0. In these regions measurement 
errors were as high as ± 10 percent of the overall AT. 

6. Improvements in surface emissivity of the heater and side 
walls may decrease measurement errors in these regions. 
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Appendix— Inframetrics Thermal Imaging System Specifications 


Explanation of Specifications 

The factors which determine the performance of the IR 
imaging system are described here. Some of this basic 
information can be found, although lacking detail and clarity, 
in the operating literature (ref. 13) of the Model 600. A 
complete list of the specifications for the Model 600 is also 
included in this appendix. 

Detector type — The Model 600 uses a single element mercury 
cadmium telluride (HgCdTe) detector which is liquid nitrogen 
cooled (77 K). The detector is mounted inside an evacuated 
metal cryogenic dewar where it is in direct contact with the 
liquid nitrogen (fig. 25), which makes the detector sensitive 
to radiation in the 8- to 14-jtm region. The response of the 
detector over this spectral bandpass is shown in figure 26. 


Thermal sensitivity .— The thermal sensitivity can be 
characterized by two parameters— the minimum resolvable 
temperature difference (MRTD) and the noise equivalent 
temperature difference (NETD). The MRTD is appropriate 
to define the thermal sensitivity of the system in terms of small 
temperature differences which are displayed on the monitor. 
The MRTD is determined by a test pattern of small vertical 
bars of different temperatures, which are seen on the monitor 
by a human observer. The temperature difference between the 
bars is decreased until the observer can no longer distinguish 
a contrast. This resulting temperature difference is the MRTD. 
The NETD is used to characterize the thermal sensitivity of 
a quantitative imager. The NETD is the temperature difference 
between two objects that results in a signal equal to the random 
background noise of the system. The MRTD and the NETD 
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(a) Scanner head. 

(b) Electronics and control unit. 

Figure 25.— Inframetrics Model 600 block diagrams (reprinted by permission of Inframetrics, Inc.). 
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Figure 26.— Detector spectral bandpass— relative response of Model 600 

system (reprinted by permission of Inframetrics, Inc.). 

for the Model 600 are 0.1 and 0.2 °C, respectively. The 
NETD can be reduced to 0.05 °C with the use of the image 
averager option. 

Image Resolution.— To describe image resolution, field of 
view (FOV), instantaneous field of view (IFOV), and slit 
response function (SRF) must be considered. 

Field of view: The FOV of the Model 600 is described in 
angular terms— 20° horizontal by 15° vertical. Therefore the 
size of the object that can be viewed is dependent on the 
distance from the scanner to the optics. 

Instantaneous field of view : The Model 600 is a scanning 
imager because a smaller field of view or instantaneous field 
of view, commonly known as the spotsize, is scanned over 
the total FOV. The radiation from each IFOV is focused on 
the detector; as the IFOV is scanned over the FOV the 
complete image is constructed. The size of the IFOV is 
determined by the characteristics of the optics and detector. 
The IFOV size (fig. 27) for the Model 600 is 2 by 2 mrad 
(0. 1 1 0 by 0. 1 1 °), or 0.08 by 0.08 cm, at a distance of 40 cm. 
An IFOV of this size allows 175 horizontal by 130 vertical 
IFOV’s per image. After being constructed in the manner 
described above, each image is digitized into pixels (256 by 
256 for the Model 600), which are smaller than the detector 
elements (IFOV’s). In this way the pixel digitization does not 
reduce the system performance. One must realize that a 256 
by 256 pixel image does not contain nearly that many true 
image resolution elements (IFOV’S) (ref. 13); therefore the 
resolution of the image is determined by the scanner, not 
digitization by the electronics. 



Figure 27.— Fields of view: FOV is 15° vertical by 20° horizontal and IFOV 
is 2 by 2 mrad. 


Slit response function : The spatial resolution is defined by 
the SRF. Figure 28 shows a three-dimensional plot of the 
response of a typical imaging system to a point radiation 
source, where the vertical axis represents intensity 
(temperature) and the x- and y-axes are spatial representations. 
Notice that the response of the system does not correspond 
to the step change in intensity (fig. 29). The system response 
to the input is less than 100 percent; that is, the reported 
intensity is less than 100 percent of the amount emitted by the 
target. This happens because the detector cannot fully respond 
as the IFOV is scanned over a point source. As this point 
source becomes larger, the response approaches 100 percent. 
It is convenient to characterize this aspect of a thermal imaging 
system by the variation in response to a vertical slit of varying 
widths, which is essentially a step change in temperature over 
the slit. The curve of the variation in response as a function 
of slit width is called the slit response function (fig. 30). The 
measured temperature of an object whose approximate SRF 
is known is a function of the object temperature and average 
surrounding temperature; that is, 

^apparent = SRF X T 0 bject 4” (1 SRF) X T surroun d 
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Figure 28.— Three-dimensional imager response to a point radiation source 
(reprinted by permission of Inframetrics, Inc.). 
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Figure 29.— Imager response to a step change in intensity. 
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Figure 30.— Slit response function with respect to slit width (reprinted by 
permission of Inframeterics, Inc.). 

The error between apparent and true temperature is greater 
when a large difference exists between the object and 
surrounding temperatures (ref. 13). The Model 600 detector 
subtends a 2 by 2 mrad solid angle, which is the IFOV. When 
the vertical slit width is equal to the detector angle (2 mrad), 
the detector output is approximately 60 percent of the actual 
temperature. (See fig. 30.) This causes some error in the 
reported temperature. Industry convention suggests that a 
target be five times the size of the IFOV in order to produce 
an accurate temperature (ref. 14). If there are sharp thermal 
gradients present on the target, some error in the reported 
temperature will result because of the response of the system. 

Information rate .— The information rate depends on image 
acquisition, image display, and the scanning mechansim. 
Image acquisition : Scan speed is the parameter that 
determines how quickly the instrument will respond to temper- 
ature changes. The scan speed is the rate at which complete 
thermal images, or full two-dimensional temperature fields, 
are updated. Each frame is acquired at intervals of 1/60 sec, 
which is called the vertical scan speed. To acquire a complete 
image the Model 600 scans 130 horizontal lines at 125 /isec 
per line or the entire field of view in 1/60 sec (60 Hz). The 
time required to scan the IFOV’s in each line of the frame 
is called the horizontal scan speed. The output of the detector 
is continuous as the IFOV is scanned over each line in the 
FOV. This continuous signal comprises the complete thermal 
image, which is then digitized by the electronics. After digital 
conversion and processing, the digital image is converted to, 
and displayed as, a television industry standard analog signal 
(EIA RS-170 or NTSC composite). 

Image display : The format of the Model 600 displayed image 
is defined by television industry standards. This format allows 
the Model 600 to output a TV signal which can be viewed 
on a conventional TV monitor. The detected image is digitized 
into 256 horizontal by 100 vertical samples. The images are 
processed and converted back to, and displayed as, analog TV 
signals (100 horizontal lines displayed at 1/60 sec). Two 
successive frames are interlaced (combined by alternating TV 
lines from each frame) to produce a picture with 200 vertical 
lines. The resulting frame rate, the rate at which each image 
is updated on the monitor, is 1/30 sec. Because of the 
interlacing, the response time of the Model 600 is 1/30 sec. 


Even though the horizontal scan speed of the displayed image, 
53 /isec per line, is faster than the horizontal scan speed of 
the acquired image, the horizontal acquisition rate of 125 /xsec 
per line is the limiting factor; thus the acquisition rate 
determines the response time of the instrument. 

Scanning mechanism : The Model 600 uses electro- 
mechanical galvanometers to perform the horizontal and 
vertical scanning. A magnetic field rotates a mirror back and 
forth several degrees against a stiff spring under closed-loop 
electrical control. The horizontal scanning is done in a 
sinusoidal mode, and the vertical scanning, in a sawtooth mode 
(ref. 13). The scanning mechanism is enclosed in an evacuated 
module that reduces acoustic noise caused by the rotating 
hardware (see fig. 26). 

Dynamic Range.— The dynamic range of the Model 600 is 
7 bits. After digitization, each pixel in the image is coded into 
128 discrete gray levels. The smallest temperature difference, 
at higher temperature ranges, that can be detected quantita- 
tively by the Model 600 is determined by the dynamic range 
and the selected temperature range. If, for example, the 
temperature range selected is 100 °C, then there would be 
100/128 °C per discrete gray level (0.8 °C sensitivity). The 
dynamic range can be extended to 8 bits (256 discrete gray 
levels) if the image averager option, which averages either 
2, 4, or 16 successive frames, is used. 

Calculation of absolute temperatures . —Because the Model 
600 only measures radiation emitted from a target, the absolute 
temperature must be calculated from this information. The 
relationship between radiation received by the Model 600 and 
the temperature of the radiating body is determined empirically 
through unique sets of calibration curves (which correlate 
received radiation with target temperature) and by viewing a 
reference source. An internal microprocessor selects the proper 
calibration curve for the current operating conditions. The 
calibration curves for each unit are determined by the 
manufacturer prior to shipment of the equipment. Before the 
conversion to absolute temperature takes place, the radiometric 
equation (ref. 14) must first be solved to account for factors 
such as target emissivity e t and background radiation. With 
the assumption that the reference source, optics, and back- 
ground temperatures are all equal, the Model 600 measures 

A/ =rJe,/, + (l -e,)/*] -7 r 

where A I is the difference between apparent radiation 
intensities of the target and the reference source. The bracketed 
quantity is the sum of the radiation intensity emitted from the 
target and the background radition intensity reflected by the 
target, and r w is the transmission of the optics. The reference 
used by the Model 600 is sampled internally at approximately 
60-Hz intervals. The target emissivity and the background 
temperature are user input. After A I is calculated, the absolute 
temperature is determined from the appropriate calibration 
curve on the basis of the selected temperature range. (A 
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graphic representation of a typical calibration curve is shown 
in fig. 31.) The process is repeated for each pixel in the 
digitized image of every frame. Therefore, the Model 600 is 
able to display the temperature of each pixel at intervals of 
1/30 sec. 
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Figure 31.— Typical calibration curve where l v l n and T r are known and T t 
is calculated. 


Inframetrics Model 600 

SPECIFICATION SUMMARY 1 

SCANNER TYPE AND SPECTRAL BANDPASS 

Scanner Type: Dual Galvanometer Reflective in sealed enclosure 

Spectral Bandpass (Model 600, 600L): 8-14 pin i (Std.) 3-6 pm (Opt.), 3-14 pirn (Opt.) 

Spectral Bandpass (Model 610) 3-6 pi m and 8-14 pirn (two channels) simultaneous 

Infrared Detector Type: HgCdTe; Liquid Nitrogen Cooling 

THERMAL SENSITIVITY 

Noise Equivalent Temperature Difference: 0.2 °C (8-14 pm), 0.4 °C (3-6 fim) 

Minimum Detectable Temperature Difference: 0.1 °C (8-14 pun), 0.2 °C (3-6 pi m) 

Integrated measurements NETD: 0.05 °C (8-14 pm), 0.1 °C (3-6 pirn) 

SCAN SPEED 

Scan Rate (Image Acquisition): 4 KHz Horizontal, 60 Hz Vertical 

Output Rate (Image Display): 15 750 Hz Horizontal, 60 Hz Vertical 


IMAGE RESOLUTION 


Total Scanned Field of View: 15° Vertical by 20° Horizontal 

Electro-Optical Zoom: 8X continuously adjustable (4X for Model 600L) 


Resolution Parameter 

600, 610 (8-14 fim) 

600L (8-14 M ni) 

600,610 (3-6 

Horizontal IFOV 

2 mRad 

3 mRad 

4 mRad 

IFOV’s per Line 

175 

116 

87 

Pixels per Line 

256 

256 

256 

Vertial IFOV 

2 mRad 

3 mRad 

4 mRad 

IFOV’s per field 

131 

87 

65 

Lines per Frame 

200 

200 

200 


fim) 


Scan Linearity: Better than 5%, typical 2% Horizontal and Vertical 

Output Format: EIA RS-170 (RS-330) Monochrome, NTSC Color, 2:1 Interlace 

(3.58 MHz color subcarrier is not locked to H sync) 


1 Reprinted with permission granted by Inframetrics, Inc., North Billerica, MA 01862-2598, from Installation and Operation of the Inframetrics Model 600 
IR Imaging Radiometer. 
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SPECIFICATION SUMMARY 


Dynamic Range: 
Calibration: 

Temperature Ranges: 

Center Temperature 
Adjustment Range: 
Temperature Readout: 
i Readout Accuracy: 


MEASUREMENT PERFORMANCE 

7 bits, 128 levels, 42 dB (x2 with Image Averager) 

Internal reference sampled 60 times per second 
Calibration updated every 1/2 second 
5, 10, 20, 50, 100, 200 (Normal Range) 

500, 1000, Max. (Extended Range) 

-20° to +400 °C (Normal Range) 

0° to > +1000 °C (Extended Range) 

3 digits 

±2% of (reading -25 °C) or 1 °C (typical) 

±5% of (reading -25 °C) or 2 °C (maximum) 

User Field Calibration process restores system to typical. Specs, 
are for standard (8-14 jim) system without filters and at ambient 
of 25 °C. Typical performance degrades to maximum at extremes 
of ambient operating temperatures. 


PHYSICAL/ELECTRICAL PROPERTIES (600, 600L) 


i 

! Power Requirements: 


A.C. Power Supplies: 
Size and Weight: 


Operating Ambient: 


System; 11-17 VDC, 10 Watts 
Monitor; 11-15 VDC, 15 Watts 
Image Averager option adds 5 Watts to system 
108-125 VAC @ less than 25 Watts each 
Scanner; 8.1x4. 9x4. 8 in., 6.5 lbs. 

Electronics; 9.3x5.4x10.3 in., 7.5 lbs 
Monitor; 12.7x4.9x6.3 in., 8.5 lbs. 

-20 to +5 °C (Model 600), 0 °C to +40 °C (Model 600L) 


PHYSICAL/ELECTRONICS PROPERTIES (610) 


Power Requirements: 


A.C. Power Supplies: 
Size and Weight: 
(Complete) 

Size and Weight 
(Discrete Units) 

Operating Ambient: 


System; 11-17 VDC, 20 Watts 

Monitors; 11-15 VDC, 35 Watts X2 

Image Averager option adds 10 Watts to system 
Total 115 VAC @ less than 200 Watts 

Scanner; 9.8x4. 9x4. 8 in., 6.5 lbs. 

Control Electronics Units; 27.8x25x25.7 in. 106 lbs. 
Electronics Unit 19.0x7.0x10.2 in., 15 lbs. 

Control Panel 19.0x7.0x5.0 in., 4 lbs. 

Monitors 19.0x8.0x15.8 in., 57 lbs. 

-20 to +50 °C 


COMPUTER INTERFACE (600, 610) 


Type: 

Speed: 

Protection 


RS-232c bidirectional Serial (Dual for 610) 
4800 BAUD 
Optically isolated input 
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